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The statistical mechanical approach to Casimir problems for di- 
electrics separated by a vacuum gap turns out to be compact and ef- 
fective. A central ingredient of this method is the effect of interacting 
fluctuating dipole moments of the polarizable particles. At arbitrary 
temperature the path integral formulation of quantized particles, de- 
veloped by H0ye-Stell and others, is needed. At high temperature - the 
limit considered in the present paper - the classical theory is however 
sufficient. Our present theory is related to an idea put forward earlier 
by Jancovici and Samaj (2004), namely to evaluate the Casimir force 
between parallel plates invoking an electronic plasma model and the 
Debye-Hiickel theory for electrolytes. Their result was recently recov- 
ered by Hpye (2008), using a related statistical mechanical method. In 
the present paper we generalize this by including a constant permittiv- 
ity in the description. The present paper generalizes our earlier theory 
for parallel plates (1998), as well as for spherical dielectrics (2001). We 
also consider the Casimir force between a polarizable particle and a 
conductor with a small density of charges, finding agreement with the 
result recently derived by Pitaevskii (2008). 
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1 Introduction 



The large interest in the Casimir effect in recent years has resulted in a flood 
of papers, most of which are concerned with the field theoretical approach 
to the problem. This is quite natural, and is in accordance with the spirit in 
the original paper of Casimir paper [TJ . Some general reviews of the Casimir 
effect can be found in Refs. [2j [3l HI [5] . Quantum field theoretical methods 
have proved to be quite effective even in the presence of dielectric media, 
although there are fundamental problems here, especially as regards perfect 
conductor boundary conditions at sharp surfaces [6]. 

Now there exists an alternative and probably less known alternative route 
to derive a theory for the Casimir effect, namely to start from statistical 
mechanics and regard the effect to be due to interacting fluctuating dipole 
moments of polarizable particles. As shown by Brevik and H0ye, the Casimir 
force between a pair of polarizable point particles can be recovered [7j. In 
this context, the path integral formulation of quantized particle systems was 
utilized. The method had earlier been applied to polarizable fluid systems, by 
H0ye, Stell, and others [8]. Subsequent generalizations led to the well known 
Lifshitz formula for parallel plates [9]. Similar evaluations were performed in 
Refs. [mig. 

We wish to emphasize, as a general remark, that the quantum statistical 
mechanical method is quite compact and effective. For instance, the Casimir 
free energy due to the dispersion force between two polarizable particles as 
given by Eq. (5.15) in [7] is derived in a way which according to our view is 
simpler than the conventional field theoretical method. The latter involves 
fourth order approximation of perturbation theory [T3]. Moreover, we ought 
to point out that the statistical mechanical approach opens new perspec- 
tives in the sense that one avoids quantizing the electromagnetic field itself. 
Instead, the field is playing the role of an agent permitting the interaction 
between polarizable particles. The two pictures are physically equivalent 

m U ED E2]. 

The topic with which we will be concerned with below, is related to earlier 
work of Jancovici and Samaj [2] . They realized that it should be possible to 
evaluate the Casimir force between metallic plates by employing an electronic 
plasma model. They considered accordingly charged particles at low density 
in a neutralized background, limiting themselves to the classical (i.e., the 
high temperature limit) for which the Debye-Huckel theory of electrolytes is 
fully applicable. Calculating the pair correlation function, and from that the 
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local ionic density at the surfaces of the plates, they were able to recover 
the conventional Casimir result for real metals in the high temperature limit. 
The issue of temperature corrections to the Casimir force has been subject 
to a lively discussion in the contemporary literature [15J. The ionic plasma 
has also been extended to the quantum mechanical case by use of the path 
integral formalism from a statistical mechanical viewpoint |16j . 

The recent paper of H0ye [T7| reconsidered the ionic plasma in the classical 
limit, using the statistical mechanical method in a different way to obtain 
the Casimir force. The correlation function was used to directly evaluate the 
average force between pairs of particles in the two plates and then the total 
force was found by integration. This is the same method as used in Refs. [HE] ■ 
The result was found to be in agreement with Ref. |14j . A notable feature 
of this approach is that it demonstrates how the modification of the density 
profile at the surface is a perturbing effect that can be neglected to leading 
order. 

In Ref. [T7] the two parallel plates (or slabs) were assumed to be filled with 
ions, but no dielectric properties due to the medium itself were envisaged. 
The present paper presents an extension and generalization of the theory in 
|17j . in the sense that the dielectric properties are included. That means, we 
introduce permittivities e in the slabs. The two permittivities are taken to 
be constant, and equal. Magnetic properties in the media are omitted. In the 
gap region between the plates, we assume a vacuum. Our method consists 
in replacing polarizations by equivalent charges. That is, the polarization P 
is taken equivalent to an induced charge density 

P/ = -V-P. (1) 

This follows from the connection between the electric field E, the polarization 
P, and the charge density p c , 

V ■ E + 4ttV • P = 47rp c . (2) 

When evaluating the Casimir force based upon the ionic fluid correlation 
function, the induced charge can be considered as a free charge. Accordingly, 
no separate charged dipole and dipole-dipole correlation functions are needed. 
In addition, the dielectric properties are expressed via the influence of the 
permittivity upon the charge-charge correlation function. 

In the next section we establish the basic formalism for the pair corre- 
lation function and thereafter, in section 3, derive the Casimir force as a 
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generalization of the expression found earlier for the case where charges were 
absent [5J- In section 4 we consider the force between a polarizable parti- 
cle and a slab (half-space) containing free charges. The treatment in that 
section is related to recent work of Pitaevskii [18]. Finally, in section 5 we 
consider the Casimir interaction between spherical dielectric shells and ex- 
tend the theory given earlier in Ref. [TT] to the case where there are free 
charges present. 

2 The correlation function 

To obtain the correlation function we start from the Ornstein-Zernike (OZ) 
equation [17] 



here extended to non-homogeneous fluids. In this equation h(v2,ri) is the 
correlation function, c(r 2 , ri) is the direct correlation function, and p(r') is 
the number density of charges. 

It may be appropriate to give a brief account of the background for this 
equation. The correlation function or pair correlation function h(r2, ri) is 
related to the pair distribution function g(r2,ri) through 



The g(r 2 , Ti) is the probability density for one particle to occupy the position 
r 2 while another particle occupies the position ri. For an ideal gas with 
uncorrelated particles, g = p(ri)p(r 2 ). Thus the h expresses deviations from 
the ideal gas value. For a uniform fluid p(r) = p =constant, and h — > 
h(r 2 — ri). The above equation as originally introduced by Ornstein and 
Zernike [19J serves as a definition of the direct correlation function c(r2,ri). 
In their investigation they noted that the direct correlation function c for a 
fluid was closely related to the pair interaction itself. In view of this, Eq. (J3]) 
can be given a simple interpretation: The resulting pair correlation function 
h is the result of a direct correlation c plus correlations via other particles as 
expressed by the integral. 

For weak long range forces, the c(r, r') is to leading order related to the 
interaction i[) in a simple way [20J , 



h(T2,T 1 )=c(T 2 ,T 1 )+ / c(r 2 ,r')p(r / )/i(r',r 1 )dr / , 



(3) 



p(ri)p(r 2 )fc(r 2 , ri) = g(r 2 , r x ) - p(ri)p(r 2 ). 




/3V(r,r') 



(4) 



4 



where (3 = l//cgT, T is the temperature, and ks is Boltzmann's constant. 
This follows from the 7-ordering studied in Ref. [20j, where 7 is the inverse 
range of interaction and the limit 7 — > is considered. In the present case this 
becomes exact for large separations r — r' with particles located in separate 
plates. However, for particles at close separation (in the same plate) there will 
be deviations. For low densities they can be neglected, in accordance with 
the Debye-Hiickel theory for electrolytes. For higher densities such deviations 
will mainly change the inverse Debye-shielding length. This contributes only 
to a minor change of the effective separation between the plates [17]. We 
find no reason to consider this further here, especially since we consider 
semiconductors with low density of charges. 

As the interaction follows from the electrostatic potential between two 
charges, we have 

V 2 c(r,r')=47r/?g c 2 -i^5(r-r'), (5) 

with s(r') the permittivity at position r', and q c the ionic charge assuming 
one component for simplicity. It is then assumed that the ions are neutralized 
by a uniform background of counter ions. 
Equation ([3]) can now be rewritten 

V 2 $ - 47r/?g c 2 -Lp(r)$ = -47r—-S(r - r ), 
e(r) e(r„) 

Mr,ro) = -/3g 2 $, (6) 

where <3> is the electrostatic potential. Here we have replaced r2 and ri by r 
and ro, respectively. With parallel plates the number density is 

( p, z<0 

p(r) = <^ 0, < z < a (7) 
y p, a < z 

with equal densities (p constant) in the two media. By Fourier transform in 
the x and y directions, Eq. (jfij) becomes 

(£-*-*)* = -4^<*-*), (8) 

where the hat denotes Fourier transform. Further, we have introduced the 
quantity k £ , defined by 

n 2 £ = 4irPq 2 cP /s (9) 
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in the two media and k £ = in the vacuum gap < z < a. Physically, k £ is 
the inverse Debye-Huckel shielding length. The symbol is the wave vector 
transverse to the z direction. The solution of Eq. (jSj) can be written in the 
form 

f i e ~ ?KZ + jBe9K2 ' z o< z <0 
$ = 2ne qKZ0 I Ce-* z + C x e qz , 0<z<a (10) 
{ De~ gKZ } a< z 

where q = k±, q K = \/k 2 L + k 2 . [ For z < Zq, the — q K {z — zq) is to be replaced 
by —q K \ z — z o I m the first line of Eq. (ilDj) .] As the boundary conditions 
require that $ and ed^/dz are continuous, one finds for the coefficient of 
main interest 

D _ Aqe^~^ a A _ ( eq K - q 

' (eq K + q) 2 (l- Ae- 2 i a )' \eq K + q 

With this the pair correlation function for free charges, for zq < and z > a, 
is 

h(k ± , z, z ) = -2tx f3q 2 c De- q ^ z - Zo) . (12) 




3 Casimir force 

To obtain the Casimir force the ionic interaction ip = ip(r) in vacuum is 
needed. This is the Coulomb potential ip = q 2 /r. Its full Fourier transform 
is V> = Anql/k 2 consistent with Eq. (jSJ) with e — 1. With k 2 = k\ + k\ this 
can be transformed backwards to obtain (recall that q = k±) 

e -q\z-zo\ 

^(k ± , z-z ) = 2nq 2 c . (13) 

1 

With q 2 = k\ = k 2 x + k 2 , dk x dk y = 2nqdq, we now obtain for the Casimir 
force per unit area, for the free ions, 

p 2 [- 

f = — h(q, z, Zo)if> g (q, z - z )qdqdzdz , (14) 

where ip' z = dip/dz, the hat denoting Fourier transform with respect to the 
x and y coordinates. (We have used J fgdxdy = J fgdk x dk y / (2ti) 2 and the 
translational symmetry along the xy plane.) Cf. Eq. (6) in Ref. [17] . Equa- 
tion (114p follows from from the general structure of the expressions for the 
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Casimir free energy and force established in Refs. [7] and |9j. In those refer- 
ences the interaction was the dipole-dipole interaction while the correlation 
function was the dipole-dipole correlation function. In the present case, by 
contrast, those quantities are replaced by the charge-charge interaction and 
the charge-charge correlation function. 

With z — z Q = v,\ + U2 + a, k 2 = en 2 we obtain 

2 roc poo poo 

f = — (-2TT(3q 2 c )D(2Trq 2 ) / / e-^ +q ^ Ul+U2+a) d Ul du 2 qdq 
2tt J Jo Jo 

K A roo De -(q K + q )a 

= ~^75 / "7 —^T qdq - 15 

8?r/3 Jo {q K + q) 2 

However, this expression is missing the contribution to the force due to fluc- 
tuating dipole moments. In the evaluation of the pair correlation function we 
avoided this problem by use of the permittivity. Here, we may use the free 
charge picture of dipole moments where the induced charge density in terms 
of polarization is given by the expression (pQ). Alternatively, we can use the 
method developed by H0ye and Stell in their analysis of the ion-dipole fluid 
mixture where all the correlations functions between ions and dipoles were 
obtained for a uniform fluid |21j . 

The polarization of a particle at position r can be written as P = Sj5(r — 
rj), where Sj is the dipolar moment. The induced charge pi as given by Eq. (pQ) 
will replace the q c factors in both the interaction ffT3l) and the correlation 
function (fl2j) . But this charge pi = — V ■ P is distributed in space and 
should be convoluted with the point particle interaction ijj and the correlation 
function h. Taking into account the Fourier transform in the transverse 
directions and the exponentials in the z direction, together with the derivative 
of the delta-function, we find that at each end of ip and h the q c should be 
replaced by 

g c ^h-Sj, h = {ik x ,ik y ,±q} 

and 

q c -> h K ■ Si, h K = {ik x , ik y , ±q K }, (16) 

respectively. [For convenience we keep the notation of Ref. [9]. Thus this h 
should not be confused with the correlation function h introduced in Eq. (j3J).] 
Now one needs thermal averages. There is no correlation between ionic 
charges and dipole moments Sj for the reference system (i.e., with the ionic 
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interaction absent). For polar or polarizable particles we assume isotropy so 
that ^ 

(O = (*iy) = (C) = 3. (SixSiy) = 0,etc, 

where the hat denotes unit vectors. To obtain the Casimir force we need the 
average 

Q = ((h-s l )(h K -s l r) = i(h-h:) 

= \{kl + qq,) = \q{qn + q) (17) 

{k 2 + k 2 = k\ = q 2 ). Note that this is the same as the square root of the 
result (6.27) of Ref. [9j, with q £ replaced by q K . Further, the dipolar density 
Pd replaces the ionic density p. We may now define the quantity y via 

Ait 

3y = y/3pdk 2 > = 4vrp d «, (18) 

where a is the polarizability. Thus for the dipoles we altogether have the 
replacement 

An(3pq 2 c -> A-K{3p d (s 2 )Q = 9yQ. (19) 

This should be modified, however, by another factor (e — l)/3y, which is 
present in Eq. (5.5) in Ref. [5] or in Eq. (33) in Ref. |21j. This may be regarded 
as a contribution from correlations between neighboring dipole moments in 
a reference system where the ideal dipole-dipole interaction (times —0) has 
been subtracted from the direct correlation function. Thus we have 

47r(3pq 2 - 3(e - 1)Q = (e - l)q(q K + q). (20) 

From the free ions we have 4ir/3pq 2 = k 2 = en 2 . Adding the contribution 
(|20p we obtain 

k 2 -> ek 2 + (e- l)q(q K + q) 

= s(q 2 K - q 2 ) + {e- l)q{q K + q) = (eq K - q){q K + q). (21) 

Altogether, the resulting Casimir force will be the one in which the k 4 in 
Eq. f|T5|) is replaced with the square of expression fT2T|) . We find 

1 f°° Ap~ 2qa 

f = -wl T^W'*' (22) 

with A given by Eq. ( fTTI) . This expression is a simple generalization of the 
situation with absence of dielectric properties [17] . or absence of charges [9]. 
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Force between a polarizable particle and a 
half-space with charges 



This is a situation earlier considered by Pitaevskii [18]. We will reconsider 
it with the formalism developed above. To simplify, the polarizable particle 
will be considered to be confined in a thin plane or layer together with other 
particles of the same kind at low density. Then the electrostatic problem 
will be as before except for the boundary condition where now one of the 
half-spaces (slabs) is removed and replaced with the thin layer of polarizable 
particles of vanishing density not influencing the electric field. Thus, instead 
of Eq. (1T01) . the solution of the electrostatic problem becomes 

i. „ f -e~ qz + Be qz , < z < a 

$ = 2tt <^ « ' (23) 

with the polarizable particle(s) located at z = 0. From the usual boundary 
conditions one finds 

D = . (24) 

eq K + q 

The Casimir force is obtained by a modification of Eq. ([151) . First, one of the 
integrations is reduced to a thin layer of width Ad at u\ = 0. Secondly, the 
modification (I2T!) is needed with dielectric media. Thus 

4TT/3pq 2 c — > (eq K - q) {q K + q). (25) 

So for the half-space 

« 2 -> >4 = - q)(q K + q), (26) 

while for the thin layer with no free charges, q K q. Further, the permittivity 
E\ of the thin layer is directly related to the polarizability a and the density 
of dipolar particles >■ 0) via 

s l - 1 = 4vrapi. (27) 

So for the thin layer 

k 2 ^ k\ = 2{e 1 - l)q 2 = Mapxq 2 . (28) 
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Using Eqs. (IMj) - (12g]) in Eq. ([IS]) we find (with Ad = J du x ) 



f 



l^Ad 



J 

Jo 



oo 



qdq 



(piAd) 



2a r eq R - q 
13 Jo eq K + q 



e~ 2qa q 3 dq. 



(29) 



Now the number of particles per unit area is p\Ad, so the force upon each of 
them is f/(p\Ad). The corresponding interaction potential V is determined 
from dV/da = — f / (piAd) . Integration of Eq. (1291) thus gives 



which is the result obtained by Pitaevskii [15] . See also the comment of 
Geyer et al., and the reply of Pitaevskii [22] . 

5 Free energy of concentric spherical bodies 

Consider an inner sphere (ball) with radius a and an outer sphere with inner 
radius b (> a) and outer radius at infinity. The spheres have permittivity 
e. Between the radii a and b there is a vacuum gap. So far this is the 
situation considered by H_0ye et al. [UJ. in that paper general expressions 
were found, for arbitrary temperature, for nonmagnetic spheres including 
both the transverse magnetic and the transverse electric modes for nonzero 
Matsubara frequencies. In the present work we want to extend this to the 
situation where also free charges are present. As in the parallel-plates situ- 
ation we have to restrict the evaluation to the electrostatic or the classical 
high temperature case. With this limitation the resulting free energy will be 
a straightforward extension of previous results. 

Employing spherical coordinates the solution for the potential Eq. ([6]), 
can be written 



where Yi m are the spherical harmonics. The radially dependent term can be 
written as 




(30) 




(31) 




r < a 
a < r < b 
b < r. 



(32) 
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With no free charges (i.e., k = 0), Eq. §6§ is the Laplace equation and the 
functions e e and s e simplify to 

e e oc e oc and s e oc s oc r'. (33) 

However, with free charges present in the media the k £ given above in Eq. (Q 
will enter Eq. ([6]). Then the solution for the radial part becomes Riccati- 
Bessel functions with imaginary argument. As $z(r) should be finite at the 
origin and zero at infinity we can write 

s £ = rji(kr) and e £ = rh[ 1 \kr), (34) 

with 

- k 2 = k\ = ATi[3q 2 c p/e (35) 
in the media. In vacuum k = in e — > 0, and we can put 

and s = r l (36) 



r i+i 



since proportionality factors are not needed here. 

One notes that Eqs. (|32|) and ([51]) are are precisely those of Ref. [UJ 
for the situation with TM (transverse magnetic) waves where Matsubara 
frequencies K = —ifruj = 2im/f3 (n integer) were used. In that case the 
imaginary values were k = u/c in vacuum and k = yfeuo/c in the media. 
With this equivalence the free energy associated with the mutual interaction 
between the inner and outer spheres can be evaluated precisely as in Ref. [UJ . 
Since the expression (I32|) for $^ does not determine its magnitude we can use 
the method in Sees. IV and V of [11] to obtain it indirectly, and thus obtain 
the eigenvalues of interest. 

For parallel plates we were able to do this directly, through Eqs. ffT6l)- ff2Tl) 
above, to obtain the force ( l22l) . By integration with respect to the separation 
a the corresponding free energy in that case is found to be 

which identify the eigenvalues X q = Ae~ 2qa . 

To obtain the corresponding eigenvalues A e ^ with the concentric spheres 
we write, similarly as in [11] , 
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The quantity D is to be determined from the continuity of and e<&[ at the 
surfaces. This gives the equations 

e ae + Bs a£ = Ce a + C\S a , 

e(e' a£ + Bs' a£ ) = Ce' a + C lS ' a , 
Ce b + Cis b = De b£ , 

Ce' b + Cis' 6 = eDe' b£ , (39) 

where subscripts a and b denote the radial positions and the prime denotes 
differentiation with respect to r. [Due to a calculational error in [TT] the 
corresponding equations (33) in that reference differ from those above in 
that e is replaced by its inverse. The error may be related to the somewhat 
involved discussion in [11] on the conditions for the transverse magnetic field 
and the related electric field for nonzero frequencies.] Solving for D one also 
needs D to determine A £ ;. As explained in Sec. IV of [TT], this can be 
obtained by considering single potential bonds between the two spheres. We 
then first remove the inner sphere by putting a = to obtain from the two 
last members of Eq. (139]) 

D = Do = C2 C, with c 2 = e 'f b - ebS 'b (40) 

ee be s b - e b£ s' b 

Then, removal of the outer sphere by putting b = oo and C\ — yields from 
the two first members of Eq. ( 1391) 

C = Coo = Cl , with d = £{eaeS '^ ~ e 'f ae) , (41) 

from which Dq = C\C2- Finally solving the full set of equations (139]) one finds 
D. This further used in Eq. (j3"81) yields 

x (ffofL ~ s' a s a£ )(ee b e' b£ - e' b e b£ ) 
[ee a s' a£ - e' a s a£ ){ee b£ s b - e b£ s b ) 

With no free charges, k £ = 0, by which e £ = e = \/r l+l and s £ = s = r l . 
Then X £ i simplifies to 

d (el + l + l)[e(l + l) + l}\b) ' 1 1 
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as found in Ref. [UJ. The resulting free energy F is now the zero frequency 
(K = 0) TM mode term of equation (40) in [TT] . 



with expression fT4"2|) inserted, where k is given by Eq. (1331) . 

6 Summary 

With the use of a statistical mechanical method the Casimir force between 
dielectric slabs containing free ions (i.e., semiconductors) has been evaluated 
in the high temperature classical limit. Further, for high temperatures the 
free energy of interaction between a polarizable particle and a semiconductor 
slab has been obtained. Agreement with earlier results of Pitaevskii is found 
[18] . Finally, in the same limit the free energy of interaction between two 
concentric semiconducting dielectric spheres separated by a vacuum gap is 
found. 




(44) 



i=i 
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